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^ C| Abstract. In particle physics, most of the classical models consider that the centre 

Qh' of mass and centre of charge of an elementary particle, are the same point. This 

presumes some particular relationship between the charge and mass distribution, a 
feature which cannot be checked experimentally. In this paper we give three different 
kinds of arguments suggesting that, if assumed different points, the centre of charge 
of an elementary spinning particle moves in a helical motion at the speed of light, and 
it thus satisfies, in general, a fourth order differential equation. If assumed a kind of 
|>y rigid body structure, it is sufficient the description of the centre of charge to describe 

also the evolution of the centre of mass and the rotation of the body. This assumption 
of a separation betwen the centre of mass and centre of charge gives a contribution to 
the spin of the system and also justifies the existence of a magnetic moment produced 
by the relative motion of the centre of charge. This corresponds to an improved model 

fvj ■ of a charged elementary particle, than the point particle case. This means that a 

Lagrangian formalism for describing elementary spinning particles has to depend, at 

**) | least, up to the acceleration of the position of the charge, to properly obtain fourth 

order dynamical equations. This result is compared with the description of a classical 



qq ' particles 

o' 



Dirac particle obtained from a general Lagrangian formalism for describing spinning 
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1. Introduction 

It is implicitely assumed that the independent degrees of freedom of any mechanical 
system satisfy, in general, second order differential equations, as it happens for the 
Newton equations of the centre of mass of the mechanical system. This is reflected 
in the general assumption that the Lagrangian of any mechanical system is an explicit 
function L(t, q iy qj) of the time t, the n degrees of freedom q^ and their time derivatives <&. 
Lagrangians depending on higher order derivatives are scarcely used, and their usefulness 
is left to some especific problems. In this work we are going to give three different kinds 
of arguments to justify the use of Lagrangians depending up to the acceleration of the 
centre of charge for describing charged elementary spinning particles. The centre of 
charge of an elementary particle thus satisfies fourth order differential equations. 

In section |2] we analyze the motion of the centre of charge of a charged rigid body, 
considered as a model of a classical elementary particle. In section [3] we suggest that the 
order of the invariant differential equation satisfied by the evolution of a point depends 
on the number of parameters of the kinematical group of the theory. Since the Galilei and 
Poincare groups are ten-parameter groups, and if the point represents the position of the 
centre of charge of an elementary particle, the corresponding differential equation must 
be of fourth order. Section H] analyzes the most general differential equation satisfied by 
a point in three dimensional space, which is of fourth order. Finally, section [5] describes 
the model of a Dirac particle obtained from a general formalism for describing classical 
spinning particles and compares its structure with the previous predicted motions. 

2. Rigid body arguments 

Let us consider that an elementary particle is described as a nonrelativistic rigid body. 
A rigid body is a mechanical system of six degrees of freedom. Three represent the 
position of a point and the other three the orientation of a body frame attached to that 
point. Usually, it is described by the location of the centre of mass, which is represented 
by the point q, and the orientation by the principal axis of inertia located around q. 
The centre of mass satisfies second order dynamical equations and moves like a point 
of mass m, the total mass of the system, under the total external force. The torques of 
the external forces produce a change in the orientation of the body. In this way a rigid 
body moves and rotates. 

If we describe the evolution of a different point r, it will follow, in general, a helical 
trajectory around the centre of mass, like the one depicted in the figure [TJ 

If an elementary particle is a charged rigid body, it is clear that we also need 
to know its electromagnetic structure, which can be reduced to the knowledge of the 
centre of charge and the different multipoles. If assumed some spherical symmetry for 
the electric field produced by the particle, we are left with the location of the centre of 
charge and no further multipoles. The position of this point will be used to determine 
there the actions of the external fields and to compute from there the fields generated by 
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Figure 1. Motion of the centre of mass q and and another arbitrary point r of a rigid 
body. 



the particle. In general, depending how the mass and charge are distributed, these two 
points will be different points as we shall assume here. Therefore, if we try to describe 
the evolution of the centre of mass, we have to determine also at any time the location 
of the centre of charge to compute the external forces. Newton's dynamical equations 
for the centre of mass will be written as 

m^f = e (E(t, r) + ^ x B(t, r)\ = F(t, r, dr/dt). (1) 

The electromagnetic force F depends, in general, on the electric and magnetic external 
fields defined at the charge position r and on the velocity of the charge dr/dt which 
appears in the magnetic term. 

The relative motion between r and q when analyzed by an observer who sees the 
centre of mass at rest like in figure [2] is a circular motion with constant velocity, in 
particular in the free case. If we define a unit vector n in the direction of the normal 
acceleration d 2 r/dt 2 of point r, 

1 d 2 r 

n ~M~d! 1, 

where R is the radius of the circular motion, then, the centre of mass position can be 
written as 

For the centre of mass observer, the angular velocity is also orthogonal to the plane 
subtended by the velocity u = dr/dt and acceleration du/dt of point r, and given by 

1 du . . 

w = —u x —-. (3) 

u 2 dt K J 
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Figure 2. Motion of the centre of charge r around the centre of mass at rest q. 



In this way, we have also solved the problem of the rotation of the charged rigid body 
by analyzing the evolution of just the centre of charge. 

Therefore, it will be simpler, from a theoretical point of view, just to describe the 
evolution of a single point, the centre of charge r, instead of the orientation of the rigid 
body and its centre of mass q, which will be in some average position of point r, and 
obtained from (TSJ) once the trajectory of r is computed. The elimination of the d 2 q/dt 2 
among equations ([T]) and (J2J) will give us, in general, a fourth order differential equation 
for the variable r. If the rigid body is no longer free, we shall assume that its stiffness 
means that the relationship between q, co and r and their derivatives given in (J2J) and 
(J3]) still remain valid, even in the presence of an external interaction. 

Therefore, if the centre of mass and centre of charge of a rigid body are different 
points, it is sufficient to describe the evolution of the centre of charge. Then, the rigid 
body is reduced to a system of only three degrees of freedom which satisfy fourth-order 
differential equations. We see that a Lagrangian depending on the acceleration of point 
r, could reproduce such equations. The following nonrelativistic free Lagrangian 
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produces the dynamical equations for the point r 
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which can be factorized in terms of a point q, defined as in 
d 2 q 



in the form 
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a free motion for the point q which can be interpreted as the centre of mass, and a 
harmonic motion of point r around q with a constant frequency uj. It is the presence 
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of the acceleration term in the Lagrangian which produces the definition of a point q, 
different from r, moving freely. This Lagrangian describes a particle with spin [1], where 
its detailed analysis can be found in this reference. 

3. Invariance arguments 

Let us consider the trajectory r(t), t G [t l5 1 2 ] followed by a point of a mechanical 
system for an arbitrary inertial observer O. Any other inertial observer O' is related 
to the previous one by a transformation of the spacetime kinematical group such that 
their relative spacetime measurements of any spacetime event are given by 

t' = T(t, r; £1, ... , g a ), r' = R(t, r; g u ... , g a ), 

where the functions T and R define the corresponding spacetime transformation of the 
kinematical group G, of parameters (g 1; . . . ,g a ), among any two observers. Then the 
description of the trajectory of that point for observer O' is obtained from 

t'{t) = T{t, r(t); ft, ... , g a ), r'(t) = R{t, r{t); <?i, . . . , g a ), Vt e [h, t 2 ]. 

If we eliminate t as a function of t' from the first equation and substitute into the second 
we shall get 

r'(t') = r'(t';g 1 ,...,g a ). (4) 

Since observer O' is arbitrary, equation (j3j) represents the complete set of trajectories 
of the point for all inertial observers. Elimination of the a group parameters among the 
function r'(t') and their time derivatives will give us the differential equation satisfied 
by all the trajectories of the point. Let us assume that the trajectory is unrestricted 
in such a way that the above group parameters are essential in the sense that no 
smaller number of them gives the same family of trajectories. This differential equation 
is invariant by construction because it is independent of the group parameters and 
therefore independent of any inertial observer. If G is either the Galilei or Poincare 
group, it is a ten-parameter group so that we have to work out in general up to the fourth 
derivative to obtain sufficient equations to eliminate the group parameters. Therefore 
the order of the differential equation is dictated by the number of parameters and the 
structure of the kinematical group. If the point r represents the position of the centre of 
charge of an elementary particle we get again that it satisfies, in general, a fourth order 
differential equation. But at the same time it is telling us that to obtain the invariant 
differential equation satisfied by the centre of charge of an elementary particle, it is 
sufficient to obtain its trajectory in an arbitrary reference frame, and to follow the 
above procedure of elimination of the group parameters. 

4. Geometrical arguments 

A continuous and differentiable curve in three-dimensional space, r(s), is a regular curve 
if at any point s has a tangent vector r(s) ^ 0. It has associated three orthogonal unit 
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vectors, t, n and b, called respectively the tangent, normal and binormal. If using the 
arc length s as the curve parameter, they satisfy the Frenet-Serret equations 

t = KTl, n = —nt + rb, b = —rn, 

where the overdot means the derivative with respect to s. The knowledge of the 
curvature n(s) and torsion t(s), together the boundary values ^(0), t(0), n(0) and 
6(0) = t(0) x n(0), completely determine the curve, because the above equations are 
integrable. In terms of the vector uj = rt + nb, known as Darboux vector, the Frenet- 
Serret equations can be rewritten as 

i — u} x t, h = uj x n, b = uj x b, 

so that, in units of arc length, Darboux vector represents the instantaneous angular 
velocity undergone by the local triad frame. 

By changing the notation r( k \s) = d k r/ds k , we have 

7-W = t, r^ 2 ' = nn, r^ 3 ' = kn + n{—Kt + rb) 

and this allows us to eliminate the three unit vectors t, n and 6, in terms of the three 
derivatives r^ k \ k = 1,2,3. If we replace them in the next order derivative r^ A \ one 
obtains the most general differential equation satisfied by the point r, i.e., the fourth 
order differential system 
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Let us consider that an elementary particle, instead of being a rigid body, is just a 
localized mechanical system. By localized we mean that, at least, it is described by the 
evolution of a single point r. This point could be the centre of mass, but, as mentioned 
before, in order to determine the external forces, we also need to know the location 
of the centre of charge to compute the actions of the external fields. Let us assume 
that the elementary particle is charged. By the previous arguments, if its electric field 
is spherically symmetric, we are reduced to know the evolution just of the centre of 
charge. The particle will have a centre of mass but we make the assumption that the 
centre of mass and the centre of charge are not necessarily the same point. 

Let us consider that the geometrical regular curve r(s) in three-dimensional space 
represents the trajectory of the centre of charge of an elementary particle. When the 
corresponding inertial observer uses its time as the evolution parameter, kinematics 
enters into the scene. Let us assume now that the motion of the particle is free. This 
means that we cannot distinguish one point of the evolution from another, so that the 
motion has to be at a constant velocity such that the arc length ds = \u\dt, where 
u = dr/dt is the velocity of the charge, must be independent of the time t. Otherwise, 
if ds is not the same we can distinguish one instant of the evolution from another, as 
far as the displacement of the charge is concerned. At the same time, Darboux vector 
has to be also independent of time. The Frenet-Serret triad moves and rotates in a 
free motion with constant linear and angular velocities. The curvature and torsion are 
necessarily constants of the motion. 
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Thus k = f = 0, and, in the free case, the equations (J5J) are reduced to 

r (4) + («2 + r 2) r (2) = _g_ ( r (2) + ^2 + r 2) ^ = Q 

If the curvature and torsion are constant the curve is a helix, which can be factorized 
in terms of a central point 

1 d 2 r d 2 q 2 2 2 

which is moving in a straight trajectory, while the point r satisfies 

£+<->-,) = «>, 

an isotropic harmonic motion of frequency u>, around point q. The point q clearly 
represents the centre of mass position of the free particle and its expression in terms of 
the centre of charge position is exactly the same as in the case for the rigid body (J2J). 

The centre of charge of a free elementary particle is describing a helix at a constant 
velocity for any inertial observer. If we make a nonrelativistic analysis, the relationship 
of the velocity measurements among two arbitrary inertial observers O and O', is given 
by u' = Ru + v, where v is the constant velocity of O as measured by O' and the 
constant rotation matrix R is their relative orientation. Now, 

u ' 2 = u 2 + v 2 + 2v Ru. (6) 

In a relativistic analysis 

2 

/ Ru + 7V + jj2^(v ■ Ru)v 



u 



7=(i-«7* , r 1/a , 



7(1 + v ■ Ru/c 2 ) 

where v is also the velocity of observer O measured by O' and R represents their relative 
orientation, and thus 



m U — C 2 



u 2 = 2 + c ■ (7) 

7 2 (1 + v ■ Ru/c 2 ) 2 

Taking the time derivative of both expressions ([6]) and (J7]) if v! has to be also constant 
for observer O', we get that v ■ Rii = 0, irrespective of v and of the rotation matrix 
R. This means that the vector u must be a constant vector. The centre of charge 
necessarily moves along a straight trajectory at a constant velocity, for every inertial 
observer, and the above general helix degenerates into a straight line and because the 
point is not accelerated q = r. 

This is the usual description of the spinless or pointlike free elementary charged 
particle, whose centre of charge and centre of mass are represented by the same point. 

However, in this relativistic analysis, there is one alternative not included in the 
nonrelativistic approach. The possibility that the charge of an elementary particle will 
be moving at the speed of light and, in that case, u — v! — c, for any inertial observer. 
This means that the centre of the helix is always moving at a velocity \dq/dt\ < c, and, 
if it represents the centre of mass, this particle is a massive particle. In a variational 
description of this system the Lagrangian should depend up to the acceleration of the 
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point r in order to obtain fourth order differential equations. This dependence on the 
acceleration will give a contribution to the spin of the particle and there is also another 
contribution from the rotation of the system, because the body frame rotates with 
angular velocity w. The motion of the charge around the centre of mass produces the 
magnetic moment of the particle. 

When we introduce the time as the curve parameter if the trajectory is necessarily a 
helix at a constant velocity u, it is equivalent the use of the arc length s = ut or the time 
as the curve parameter. But if the curve is regular it must have at any point a tangent 
vector r^ 1 ' 7^ 0. This implies that the velocity of the charge can never be reached by 
any inertial observer. Otherwise, the velocity of the charge will be zero at a certain time 
but different from zero for subsequent times, because the point is accelerated, which is 
contradictory with the requirement that the point describes a trajectory with a constant 
velocity. 

This kinematical argument requires that the kinematical group of spacetime 
symmetries must contain a limit velocity which cannot be reached by any inertial 
observer. Among the possible ten-dimensional kinematical groups found by Bacry and 
Levy-Leblond [2], only the De Sitter groups SO(3,2) and 50(4,1) and the Poincare 
group V contain such a limit velocity c. If the spacetime is flat then only the Poincare 
kinematical description is singled out by this requirement, so that the charge, necessarily, 
must be moving at this limit velocity. 

In summary, there are only two possibilities for a free motion of the centre of charge 
of an elementary particle. One, the charge is moving along a straight line at any constant 
velocity, the centre of mass is the same point, the system has no magnetic moment and 
it can be described either in a relativistic or a non-relativistic framework. In the other, 
the particle has spin and magnetic moment, and the charge moves along a helix at the 
speed of light. Because all known elementary particles, quarks and leptons, are spin 1/2 
particles, and the charged ones have magnetic moment, we are left only with this last 
possibility if we want to give a more improved description of an elementary particle. 
Here, only the relativistic description is allowed. 

This is consistent with Dirac's theory of the electron, because the eigenvalues of the 
components of Dirac's velocity operator are ±c. This means that Dirac's spinor ip(t,r) 
is expressed in terms of the position of the charge r, because the external fields A^it, r) 
are defined and computed at this point. 

This last possiblity is the description of the centre of charge of a relativistic spinning 
elementary particle obtained in the kinematical formalism [1], and which satisfies Dirac's 
equation when quantized. The classical structure of this Dirac particle is analyzed in 
the next section. 

In this formalism Dirac particles are localized and also orientable mechanical 
systems. By orientable we mean that we have to attach to the above point r, a local 
cartesian frame to describe its spatial orientation. This frame could be the Frenet- 
Serret triad. The rotation of the frame will also contribute to the total spin of the 
particle. When quantizing the system, the spin 1/2 is coming from the presence of 
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the orientation variables. Otherwise, if there are no orientation variables, no spin 1/2 
structure is described when quantizing the system. This twofold structure of the classical 
spin has produced a pure kinematical interpretation of the gyromagnetic ratio [3] . The 
dependence of the Lagrangian on the acceleration is necessary for the particle to have 
magnetic moment and for the separation between the centre of mass and centre of 
charge. 

5. Kinematical description of a Dirac particle 

In the kinematical formalism [1], an elementary particle is, by definition, a mechanical 
system which in addition to being indivisible, as a consequence of the atomic hypothesis 
[I], it can never be deformed so that all allowed states are only kinematical modifications 
of any one of them. This means that when the state of an elementary particle changes 
it is possible to find an inertial observer who measures the particle in the same state as 
before. An electron, if not annihilated, always remains an electron under any external 
force. This means that in a variational approach the initial x\ and final X2 states of the 
evolution are related by a transformation of the kinematical group x<i = gx\. Therefore, 
the boundary variables of the variational approach, necessarily span a homogeneous 
space of the kinematical group of space-time symmetries. When quantizing all classical 
systems characterized by such homogeneous spaces, their Hilbert space of pure states 
carries a projective unitary irreducible representation of the kinematical group [5j. It 
thus satisfies Wigner's definition of a quantum elementary particle. 

In this way, the parameters of the kinematical group become the classical variables 
we need to consider, as the boundary values of the variational formalism for describing 
an elementary particle. In the relativistic and non-relativistic approach, these variables 
are reduced to the ten variables t, r, u and a, interpreted respectively as the time, 
position of the charge, velocity of the charge and orientation. In the relativistic case 
we have three disjoint, maximal homogeneous spaces of the Poincare group spanned by 
these variables with the constraint either u < c, u = c or u > c. It is the manifold with 
u = c, as suggested by the kinematical arguments of the previous section, which leads to 
Dirac's equation when quantizing the system. Because the Lagrangian depends on the 
next order derivative of the boundary variables, it thus depends also on the acceleration 
of the point r and on the angular velocity. It is thus clear that the point r cannot be the 
centre of mass because satisfies fourth-order differential equations. Because the external 
interaction is defined at this point r, is why we consider it represents the position of the 
charge. 

Then, for a Dirac particle, the charge located at point r, is moving at the speed of 
light u = c. The classical expression which gives rise to Dirac's equation is 

7r n 1 „ fdu 

H = P-U+—S ■ l — xu 
c 2 \ at 

where the energy H is expressed as the sum of two terms, P ■ u, or translational energy 
and the other, which depends on the spin of the system, or rotational energy. The spin 
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comes from the dependence of the Lagrangian L of both, the acceleration ii, and the 
angular velocity w, and if we define 

dii ' duj ' 

it takes the form 

S = ux U + W = Z + W. 

The first part Z = u x U, or zitterbewegung part, is related to the separation between 
the centre of charge from the centre of mass and takes into account this relative orbital 
motion. It quantizes with integer values. The second part W is the rotational part 
of the body frame and quantizes with both integer and half-integer values. The total 
angular momentum with respect to the origin of observer's frame is 

J = rxP + S, 

so that the spin S is the angular momentum of the system with respect to the centre 
of charge r, and not with respect to the centre of mass q. By this reason, it is not a 
conserved quantity for a free particle, but satisfies the dynamical equation 

dS „ , . 

This is exactly the same dynamical equation satisfied by Dirac's spin operator in the 
quantum case. This has to be taken into account when comparing the analysis of this 
spin with other approaches, for instance, with Bargmann-Michel-Telegdi spin observable 
[6], which clearly represents the angular momentum with respect to the centre of mass 
of the system. Once a mechanical system has two distinguished points, the centre of 
charge and centre of mass, we must clarify with respect to which of these points is 
defined the angular momentum of the system. This is important, for instance, in the 
so called proton spin crisis. If the spin of the proton is the angular momentum with 
respect to its centre of mass, and we add the three Dirac spin operators of the three 
quarks we cannot obtain the spin of the proton. What we have to add are the three 
spins of the quarks with respect to their corresponding centre of mass, if assumed that 
the motion of the quarks is in an / = orbital angular momentum state. 

When expressed Dirac's spin and the centre of mass position in terms of the velocity 
and acceleration of the charge they take, respectively, the form 

H-u-P\du _ c 2 (H-u-P\du 

{du/dtf J ~dt XUj Q ~ r+ H\ (du/dt) 2 J ~dl 
Dirac's spin is always orthogonal to the osculator plane of the trajectory of the charge 
r, in the direction opposite to the binormal for a positive energy particle, and in the 
opposite direction for the antiparticle. This implies a difference in chirality between 
matter and antimatter. The acceleration of the charge is pointing from r to the centre 
of mass, as it corresponds to a helix. It is shown that the dynamical equation of point 
r for the free particle and in the centre of mass frame is given by 

S x u, (9) 



mc 2 
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r=Sxu/mc* 
S=Z+W 




Figure 3. Motion of the charge of the electron in the centre of mass frame. The 
magnetic moment of the particle is produced by the motion of the charge. The total 
spin S is half the value of the zitterbewegung part Z when quantizing the sytem, 
so that when expressing the magnetic moment in terms of the total spin we get a 
g = 2 gyromagnetic ratio [3]. The body frame attached to the point r, which could be 
Frenet-Serret triad, rotates with angular velocity w, has not been depicted. 



and where the spin vector S is constant in this frame, as depicted in Fig. [31 The radius 
of the zitterbewegung motion is R = S/mc, and the angular velocity u = mc 2 /S. When 
considered in the centre-of mass frame it is a system of three degrees of freedom; two 
are the x and y components of the position of the charge on the zitterbewegung plane 
and the third is the phase of the rotation of the body frame. This phase is the same as 
the phase of the orbital motion and because the velocity u = c is constant, we are just 
left with a single and independent degree of freedom, for instance, the x coordinate. 
The Dirac particle, when considered in the centre of mass frame, is equivalent to a one- 
dimensional harmonic oscillator of frequency u. If we quantize this harmonic oscillator, 
because the centre of mass is at rest the oscillator is in its ground state of energy 
hjj/2 = mc 2 . No further excited states are allowed according to the atomic principle 
|lj. The classical spin parameter S thus becomes h/2, when quantized. 

When seen from an arbitrary observer (see Figure HJ), the motion of the charge is a 
helix, so that according to (E]) Dirac's spin precess around the direction of the conserved 
linear momentum P. The spin with respect to the centre of mass is defined as 

S CM = S+(r-q)xP. 

It is a conserved quantity for a free particle. The centre of mass velocity is v = dq/dt, 
and the linear momentum is written as usual as P = ~f(v)mv. This means that the 
transversal motion of the charge is at the velocity v 'c 2 — v 2 . A moving electron takes 
a time ~f(v) times longer than for an electron at rest to complete a turn, as a result of 
the time dilation measurement. The faster the centre of mass of the electron moves the 
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Figure 4. Precession of Dirac's spin S along the linear momentum P. The tranversal 
motion of the charge takes a time "f(v) longer than when the centre of mass at rest, to 
complete a turn. The three vectors it, du/dt and —S, properly normalized, form the 
Frenet-Serret triad of the motion of the charge. The spin with respect to the centre of 
mass Scm> is a constant of the motion for the free particle. 



slower is the rotation frequency of the centre of charge around the centre of mass. The 
internal clock of a fast electron is running slower. 
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